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ABSTRACT 

Solar flare hard X-ray spectra from RHESSI arc normally interpreted in terms of purely collisional electron beam 
propagation, ignoring spatial evolution and collective effects. In this paper we present self-consistent numerical sim- 
ulations of the spatial and temporal evolution of an electron beam subject to collisional transport and beam-driven 
Langmuir wave turbulence. These wave-particle interactions represent the background plasma's response to the elec- 
tron beam propagating from the corona to chromosphere and occur on a far faster timescale than coulomb collisions. 
From these simulations we derive the mean electron flux spectrum, comparable to such spectra recovered from high 
resolution hard X-rays observations of solar flares with RHESSI. We find that a negative spectral index (i.e. a spectrum 
that increases with energy), or local minima when including the expected thermal spectral component at low energies, 
occurs in the standard thick-target model, when coulomb collisions are only considered. The inclusion of wave-particle 
interactions does not produce a local minimum, maintaining a positive spectral index. These simulations are a step 
towards a more complete treatment of electron transport in solar flares and suggest that a flat spectrum (spectral 
index of to 1) down to thermal energies maybe a better approximation instead of a sharp cut-off in the injected 
electron spectrum. 

Subject headings: Sun: flares - Sun: X-rays, gamma rays - Sun: activity -Sun: particle emission 



1. INTRODUCTION 

Hard X-ray emission has long been used as the prime 
diagnostic tool to study particle acceleration and en- 
ergy release in solar flares. From these X-ray obser- 
vations the mean electron flux spectrum (e.g. aver- 
aged over the X-ray emitting volume, see iBrown et al.l 
(|200l for de tails) can be determi ned either through for- 
ward fitting dHolman et al.ll200l or more advanced in- 
version techniques (jPiana et al Tl2fM[Ko7itar et al.ll200i 
IBrown et all 12006( 1. At higher energies, typically above 
10-20 keV, the observed hard X-ray spectrum is consid- 
ered to be due to an accelerated population of electrons 
being sto pped by the d ense chromosphere via Coulomb 
collisions (|Brownlll97ll ). The spectrum below 10-20 keV 
normally originates from the rmal coronal sources with 
temperatures of 10s MK (e.g. iKrucker fc Linll200l l. 

The Reuven Ramaty High Energy Solar Spectrome- 
ter (RHESS I) provides high resolution HXR spectra of 
solar flares (|Lin et al.ll2002f). greatly im proving on pre- 
vious measurements (Johns fc Linl 11992} ). This high en- 
ergy resolution spectra has allowed, for the first time, to 
scrutinize the X-ray and electron spectra in search for a 
non-powerlaw features, revealing vital clues about elec- 
tron acceleration and transport. In some RHESSI flares, 
the recovered mean electron flux spectrum demonstrates 
a local minima or dip between the non-thermal and 
the thermal components instead of a smooth transition 
(iPiana et al.ll2003l: iHolman et al.ll2q03l: iKasparova et all 
120051: iSui et all l2007t iKontar et all l2008al ). The pres- 
ence of the dip, as a real physical feature of the electron 
spectra, has been questioned as in many cases it can be 
attributed to phot ospheric albedo, e.g. Compton back- 
scatte red X-rays feasparo va et al.l 120051 : IKontar et al.l 
l2008a| ). However a few events hav e been found in whic h 
after isotropic albedo correction (jKontar et al.l l2008al ) , 



the X-ray spectrum is still relatively flat, so they could 
be fitted with a thick-target mo del single power -law spec- 
trum with a low energy cutoff (|Sui et al.ll2007f ) . In these 
flares a dip was not directly observed in the mean elec- 
tron spectrum, but instead inferred from forward fitting 
a model with low energy cutoff to the X-ray spectrum. 
This model has a thermal component at low energies and 
at higher energies purely collisional thick-target model of 
a single power-law of acceler ated electrons above a cutoff . 
In this thick-target scenario (jBrownll 19711 : lHolmanll2003f l . 
the dip in the mean electron spectrum originates from a 
positive slope at low energies developing below the cut- 
off as the accelerated electrons propagate from a coronal 
acceleration site downwards to the chromosphere, having 
Coulomb collisions with the background plasma. If the 
dip is real it provides important insights into flare ener- 
getics since the energy in accelerated electrons is strongly 
dependent on the low energy cutoff. 

For any reasonable X-ray producing flare, non- 
collisional beam-plasma interac tion is much faster than 
that via Coulomb collisions (jZheleznyakov fc Zaitsevl 
Il970t iKarlickvl 120091 ). Such processes are inferred to oc- 
cur in downward propagating electron beams from radio 
observations of reverse slope drift burst in flares (e.g. 
iKlein et al.lll997t lAschwanden fc Benzlll997D . Although 
generation and escape of electromagnetic radiation from 
Langmuir waves in a flaring plasma is not well under- 
stood. The role of wave-particle interactions in solar 
flares assuming stationary, time-independent injection of 
electr ons have been considered analytically and numeri- 
cally (lEmslie fc Smithlll984l:lHamilton fc Petrosianlll987l : 
lMcClementslll987D . lEmslie fe Smith! (|1984f ) have argued 
that the conditionally created distribution should be con- 
stantly being flattened by qu asi- lin ear relaxation , while 
lHamilton fc Petrosianl (|1987h and iMcClementsI (|1987l ) 
suggest that although the wave-particle interactions have 
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Fig. 1. — The background plasma density n(x) as a function of 
the height above the photosphere. Also shown is the corresponding 
plasma frequency for each of these densities. 



an important effect, the change of the electron spectra 
under stationary conditions should be minor. However, 
in a more realistic flare conditions the injection (accel- 
eration) of electrons is lik ely to be highly intermittent 
(|Tsiklauri fc Ha ruki 2008|) with a number of short du- 



ration pulses is often observed (iKiplinger et al.l 119841 ; 
iFleishman et all 119941 ; lAschwanden et al.lll998l ). so the 



time-dependent solution of particle transport equations 
accounting for wave-particle interactions should be con- 
sidered. Additionally, the previous studies did not con- 
sider the spatial and temporal evolution of the beam from 
the coronal source down into chromosphere - a crucial 
aspect when considering the propagation of an electron 
beam in comparison to X-ray observations. 

The quasi-linear relaxation process of Langmuir waves 
has been considered in higher velocity dimensions (other 
than the parallel comp onent considered here) (e.g. 
iChuraev fc Agapovlll980f) but it has only been recently 
that the 2D system has been fully numerically solved 
(jZiebell et al.ll2008f l. Even then the evolution was con- 
sidered in a spatially independent manner. In these stud- 
ies it was found that the parallel component (ID) is 
the fastest processes and likely to dominate the electron 
transport. 

In this letter we take a step towards a more com- 
plete treatment of electron transport in solar flares by 
including the spatial evolution of beam-driven Lang- 
muir wave turbulence. We numerically study the sys- 
tem self-consistcntly, simulating the propagation of an 
electron beam from the coronal acceleration site down 
to the chromosphere, considering the truncated power- 
law spectrum frequently used f or data interpretation (e.g. 
lHolmanl[20(33l : iSui et al.ll2007h . to investigate the evolu- 
tion of the mean electron flux spectrum below this cut- 
off. We demonstrate that the positive slope of the mean 
electron flux is not present, when the response of the 
background plasma via Langmuir waves to the propa- 
gating electron beam is taken into account. We also 
show that the injected electron spectrum flattens to a 
decreasing distribution due to coll ective interaction wit h 
plasma even for weak flares (e.g. lHannah et al.ll2008aD . 
Furthermore, we suggest that a flat spectrum (with spec- 
tral index to 1) down to thermal energies maybe a bet- 
ter approximation as opposed to the sharp cut-off in the 



injected electron spectrum. 

2. PARTICLE TRANSPORT AND WAVE-PARTICLE 
INTERACTION 

Following the standard model approach for interpret- 
ing solar flare hard X-ray spectra, we assume the electron 
flux spectrum of injected (flare accelerated) electrons is 
a power law, F(E) ~ E- 5 , [electrons cm s keV 1 
down to some energy Eq, typically 10-20 kcV. The initial 
ID electron distribution function (accelerated electron 
population) subsequently is also a power-law in velocity 
f(v) = F(E)/m, above a low energy cutoff vc with spec- 
tral index a = 25. For our simulations we consider such 
an initial electron distribution which is also spatially a 
gaussian of characteristic size d 
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normalised by the beam density n . 

To self-consistently follow the temporal and spatial 
evolution of an electron beam from a coronal accel- 
eration site, including the response of the thermal 
background plasma in the form of Langmuir waves, 
we use the ID equations of quasi-linear relaxation 
(IVedenov fc Velikhovl [1961 iDrummond fc Pinel [1961 
lRvutovlll969HHamilton fc Petrosianll987HKontarl2001ah 
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where f(v, x, t) is the electron distribution function [elec- 
trons cm -4 s], W(v,x,t) is the spectral energy den- 
sity [erg cm -2 ], k is the wave number of a Langmuir 
wave, n is the background plasma density and w 2 = 
4irne 2 /m is the local plasma frequency. The first com- 
ponent on the righthand side of equations @ and 
are the quasilinear terms that describes the resonant 
interaction between the electrons and Langmuir waves, 
uj p = kv. Also included arc the Coulomb collision damp- 
ening rate for both the electrons 7c F = 47re 4 nlnA/m 2 
(lEmslid H97cl and waves 7c w = 7re 4 nln A/(m 2 i4) 
(jMelrosd fl980h . Where In A = In (8 x 10 6 n -1 / 2 T) is 
the Coulomb logarithm, T is the temperature of the 
background plasma and «t = \Jk&T jm is the veloc- 
ity of a thermal electron, fee is the Boltzmann con- 
stant. Also included in equation ([3]) is the Landau 

dampening rate 7l = \/t t/ 8uj p (v/vt) exp (— v 2 /2v 2 ? ) 
(jLifshitz fc Pitaevskiilll98ll ) and the spontaneous emis- 
sion 5* = gj^ mv In (v /v T )/(4irn) (|Melrosd U980t 
iTsvtovich fc Terhaarj[l995l : lHamilton fc Petrosianlll987f ) . 

In the simulations, we take a = 7, or 5 = 3.5 from a 
cutoff of v c = 7.26 x 10 9 cm s" 1 , or E c = 15 kcV, up 



Our simulation 

9 



to maximum of vq = 2.4 x 10 10 cm s 1 
extends in velocity space from v = 7vt 
s -1 (taking T = 1MK) to v = 2.5 x 10 10 cm s" 1 . The 
initial spatial scale of the beam is d = 2 x 10 8 cm with 
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Fig. 2. — The evolution of the electron distribution f(v,x,t) (with time increasing from left to right) for the simul ation only following 
the progression of an electron beam subject to Coulomb collisions, i.e. the standard thick-target model (Brown 1971). (An animation of 
Figure [2] and [4] is available in the online journal.) 



density of n = 10 6 cm _3 . The initial number of elec- 
trons in this simulation is N as 23nod 3 = f 32 electrons, 
which is an approximation as we only have one spatial 
dimension to estimate the volume from, taking this as 
the FWHM = 2\/21n 2d, which i s typically measured 
is X-ray imaging ob servations fe.g. iHurford et all 120021 : 
IKontar et alJl2008bf ). Here we have used a modest num- 
ber of elect rons, similar to that found in a small flare, or 
microflare ([Hannah et al . 2008a. b), a typical A or B-class 
GOES flare. Since the rate of wave-particle interactions 
is proportional to electron beam density, the effects of 
wave-particle interaction will be present to a far greater 
extent in larger flares. 

We approximate the background plasma density n as- 
suming a constant of 10 10 cm -3 at coronal heights, with 
a sharp density increase at the chromosphere level, with 
further steady hydrostatic incr ease towards the photo- 
sphere ([Aschwanden et a l. 2002), see Figure[T]The initial 
electron beam is spatially centred at ho = 4 X f0 9 cm, 
see Figures [5] and S] for details. 

The equations ([2]) and ([3]) are solved numeric ally us- 
ing a f inite difference method as described by IKontar! 
(|2001bh . This is over a grid of 60 points in velocity space 
and 160 in position space. The fastest process here is 
the quasilinear relaxation, occurring on a timescale of 
r Q ~ non/u>p r~j 2 x 1Q y/nlriQ seconds. Therefore wc 
numerically solve equations {2} and ([3]) using a time-step 
at least an order of magnitude smaller. The initial spec- 
tral energy density is taken to be the thermal background 
which has reached a steady-state through Coulomb colli- 
sions and wave-particle emission/absorption. These sim- 
ulations are ran for 1 second in simulation time, enough 
time for all of the electrons to reach the highest density 
region, lose energy and then leave the simulation grid, 
joining the thermal electrons. 

2.1. Beam Coulomb Collisions 

We start by simulating the propagation of the elec- 
tron beam in the absence of waves, with only Coulomb 
collisions acting on the electrons , following the standard 
thick-target model (|Browniri971[) . Namely, only solving 
equation ([2]) and ignoring the first term on the right- 
hand side, the wave-particle interactions. The rcsult- 



10 51 

„ 10 50 
> 

y 1049 

| 10- 

LL 
> 

* 10 47 
10 46 



10 5 ' 

„ 10 50 
> 

CD 

v„ 10 49 



Beam Only 



10MK, 10 45 cm 



45 „ m -3 - 



10 

Electron Energy, E [keV] 



100 



10 4 ' 

10" 
10* 



Beam & Waves 

10MK, 10 45 cm 3 



10 

Electron Energy, E [keV] 



100 



Fig. 3. — The mean electron flux spectrum < nVF > for the 
simulation of an electron beam with only collisional dampening 
(top) and for the electron beam with the generation of Langmuir 
waves (bottom). These spatially integrated spectra are averaged 
over the lsec duration of the simulations. The dashed lines are 



a thermal model with T = 10 MK and EM 
dotted lines are the total spectrum. 



10 4 



The 



ing electron distribution f(v,x) for various times dur- 
ing the simulation is shown in Figure [5] The electrons 
with the highest velocities move quickly to lower heights 
where they encounter the sharply increasing background 
plasma density (chromosphere) below about 3 Mm. Here 
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Fig. 4. — The electron distribution f(v, x, t), top row, and wave spectral energy density W(v, x, t), bottom row, with time increasing from 
left panel to right for the simulation following the progression of an electron beam including Langmuir wave response from the background 
plasma. W is shown against electron velocity v instead of wave number k since k = ui p /v and to allow easy comparison between the two 
rows. (An animation of Figure [2] and [4] is available in the online journal.) 



the Coulomb collisions quickly cause the electrons to lose 
energy and eventually have velocities outside of the sim- 
ulation grid. At lower energies the sharp initial low 
energy cutoff is smoothed out through Coulomb colli- 
sions reducing the electrons velocity. The time aver- 
aged mean electron flux spectrum of purely collisional 
transport is shown in Figure [3J This spatially integrated 
mean electron flux spectrum < nVF > or nVF is re- 
lated to the simulated electron distribution f(v,x,t) as 
JiVF{E,t) = AJ2 [n(x)f(v,x,t)/m e ]dx where A is the 
cross-sectional area of the beam. We take this to be 
FWHM 2 m 5.5c? 2 given our ID simulation. The positive 
slope at low energies is clearly visible with the expected 
decreasing power-law above roughly the original low en- 
ergy cutoff (Figure [3] left). Ovcrplotted is an example 
model thermal spectrum to demonstrate how the local 
minima (dip) appears in the spectrum. A typical flare 
temperature of 10 MK is used, with a modest emission 
measure of 10 45 cm~ 3 to match the observations of the 
small flare we have simulated. This is a standard thick- 
target spectrum, which is used to fit and interpret X-ray 
spectrum. 

2.2. Beam-Driven Langmuir Waves 



We now follow the beam propagating with Langmuir 
waves generated by the background plasma in response 
to the beam, numerically solving both equations ([2]) and 
((3|) . The resulting electron and spectral wave density dis- 
tributions are shown in Figure [4] as a function of v and 
x for various times during the simulation. The electron 
distribution function f(v,x,t) quickly flattens to form 
a plateau-like distribution expe cted for beam-plasma in- 
terac tion via plasma waves (e.g. fZhclcznv akov fc Zaitsevl 
119701 ) . The electrons together with the waves move down 
and eventually end up in the dense regions of the atmo- 
sphere, where the transport becomes dominated by col- 
lisions. The motion of the plasma turbulence (Figured]) 
is not due to the group motion of waves, which is negli- 
gibly small dtu/dk = 3v^ e /v -C v, but appears because 
the Langmuir waves are locally generated and efficiently 
reabsorbed by the beam itself or collisionally by the sur- 
rounding plasma. The wave-particle interactions clearly 
changes the overall shape of spatially integrated electron 
flux spectrum, as shown in the right panel of Figure [3j 
Crucially, no positive slope is created in the non-thermal 
spectrum below the initial low energy cutoff, resulting in 
no dip in the overall model spectrum. This can be fur- 
ther seen when we consider the spectral index 5 of the 
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Fig. 5.— The spectral index 8(E) = -dlog(< nVF >)/dlog(E) 
of the mean electron spectrum for the simulation with the electron 
beam only subject to Coulomb collisions (black line) and with the 
beam-driven Langmuir waves (blue line) . The dashed line indicates 
the total model spectra (including thermal component) demon- 
strating a local minima in the beam only case. 

mean electron flux spectrum as a function of energy in 
Figure [S] Where as the beam only simulation produces 
a brief energy range where 5 < 0, it is always positive in 
the beam and wave case. The slightly reduced level of 
the mean electron flux spectrum also suggests that the 
generation of Langmuir waves leads to additional energy 
losses by the beam (into heating the background plasma) 
and higher number of energetic electrons will be required 
to explain the same X-ray spectra. 

3. DISCUSSION & CONCLUSIONS 

Considering both the temporal and spatial evolution, 
the mean electron flux spectrum is very sensitive to the 
generation of waves. The influence of wave-particle inter- 
actions is seen to flatten the spectral index of the elec- 
tron spectrum 5(E) = d\og(nVF)/d\og(E) below the 
break but also up to where the beam-plasma interaction 
time is faster than the electron cloud time size. The 
artificially introduced low-energy cutoff in the injected 
electron spectrum disappears not only in the local elec- 
trons distribution function but also in the spatially inte- 
grated electron spectrum once wave-particle interactions 
is taken into account. The character of beam propaga- 
tion is close to the simulation results of beam transport 
along open field lines, where collisions are normally ig- 
nored dTakakura fc Shibahashi f 976: Magolss cn fc Smith! 
[19771 IMel'Nik et al.lll999tlLi et al.ll2008fl . However, there 
are noticeable differences. In the simplistic treatment of 
spatially unifor m beam, when all but quas ilinear terms 
are ignored (e.g. iVedenov fc Velikhov|[l963l ). the genera- 
tion of waves leads to an exact plateau distribution with 
6(E) ss 0. Our simulations show that the spectral index, 
5(E) is m ore than zero, which i s the result of collisions. 
Similarly, iKontar fc Reidl (|2009l ) show that the spatially 
integrated spectrum of particles will not deviate from a 
initial power-law, but only when processes leading to ab- 
sorption of waves or removal of waves out of resonance 
are included. 



The simul ations in this letter are for typical microflarc 
parameters (jHannah et alJ l2008a| ) and given that the 
wave emission scales with the electron number den- 
sity we would expect wave-particle interactions to have 
more significant, effect in large flares. Large flares 
could constitute multiple intermittent bursts of accel- 
erated electrons directed along possibly different mag- 
netic field lines . The fast time var iations in hard X- 
ray lightcurves (jKiplinger et al.lfl984f ) indirectly support 
this idea. Numerical simulations o f reconnection sug- 
gest "bursty" electron acceleration (|Tsiklauri fc Harukil 
I2008D and spatially fra gmented electron acceleration 
( Bian fc Browning] 120081 ) which could result in electron 
propagation along the different lines. The footpoin t mo- 
tion often seen in solar flares (jKrucker et al.ll2003f ) also 
suggest that the electrons are consecutively injected onto 
field lines. In this scenario an ensemble of our simu- 
lations, multiple micro-beam injections, would lead to 
beam densities comparable to a large flare. 

The conve rgence of the magnet ic field at chromospheric 
heights (e.g. lKontar et al . 2008c) has been ignored in this 
work, however in our simulations the overall evolution of 
the energetic particles in the top part of a relatively dense 
loop, 10 10 cm -3 , is dominated by wave-particle interac- 
tions where the field is not converging. It is only in the 
denser chromosphere, where the field is likely to converge, 
that the collisions become dominant. The very fast flat- 
tening of the powerlaw distribution's low energy cutoff by 
the wave-particle interactions suggests that it is unlikely 
that such a cutoff could develop and is therefore an un- 
wise initial distribution for any model of coronally accel- 
erated electrons. The non-thermal distribution flattening 
at low energies as it transitions into the thermal distri- 
bution seems to be a realistic model. Given how strongly 
the total energy in the accelerated elections depends on 
a cutoff, or its behavi our at low energies approaching th e 
thermal distribution (|Emslidl2003tlGallowav et al.ll2005f) . 
this transition needs further study. 

The work presented here is a step towards a more com- 
plete treatment of electron transport in solar flares, mov- 
ing from the standard thick target model (ID velocity 
with collisions only) to beam-driven Langmuir turbu- 
lence (ID velocity, ID space, collisions and wave-particle 
interactions). In future work we need to extend this to 
consider 2 or higher dimensions of velocity space, the 
changing magnetic field, electron scattering as well as 
the transition between the accelerated beam and ther- 
mal distribution. 
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